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Abstract. The present paper is in a sense a continuation of [llj . 
it relies on the notation and some results. The problem tackled 
in both papers is the nature of the continued fraction expansion 
of are the partial quotients bounded or not. Numerical ex- 
periments suggest an even stronger result on the lines of Kuzmin 
statistics. Here we apply different sets of bases for the vector space 
F, where the adjunction ring 7L\\Jl!\ lives. And as a result we get 
a criterion for continued fraction convergents in terms of the coef- 
ficient vector from a lattice. 
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1. Introduction 

Stability of invariant circles in K.A.M. theory depends on the respec- 
tive rotation number. The most stable circle has = the golden 
mean ratio, as its rotation number, all partial quotients equal hi = 1. 
As in |H] boundness of partial quotients would be a bonus in represen- 
tation on a computer. Here however, experiments strongly suggest the 
opposite. Also cubic irrationals are interesting in studying quasiperi- 
odic motion [3], [9]. Here we investigate \/2 and its adjunction ring. It 
is a common belief that the partial quotients of are not bounded, 
supported by extensive computations, but no proof. 

Even more, computations suggest that their relative frequencies in the 
limit obey the Kuzmin law P{bn = k) = log2 In [12] several 

algebraic numbers were used in computations, among them ^ 
and good accordance was found with Kuzmin's statistics, for \/2 even 
too good. So later in [6] and [2J larger samples were taken and the 
anomaly seemed to disappear. We played the same game, only we had 
the advantage of more sophistical computation tools that evolved in the 
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years in between. The experimental results supporting of the stronger 
hypothesis instigated us to try towards some theoretical results. 

However questions about "big" partial quotients still linger and tease 
us. Though some explanation was given in a special case, using elliptic 
modular functions [13], [1]. 

In this paper we work towards the proof of 

Hypothesis. The partial quotients of \pl are not bounded. 

denoted as Hypothesis B in [11]. Only a very limited partial result 
is given, helping to recognize possible convergents and estimating the 
next partial quotient. 

2. Adjunction ring and ambient vector space V 

For convenience of the reader we repeat some definitions and notations 
from [llj. 

In the adjunction ring we have the unit p = 1 + \pl + and its 
inverse a = — 1 + v^, per = 1. We span the 3-dim space on p, 1, cr. 
Obviously the continued fraction expansions for \pl and a differ only 
in the starting partial quotient, so we rather consider approximations 
to cr. And in order to find ever better approximations, we construct 
series of vector space bases for V . 

Multiplicative norm is defined in TL\\pi\. Let x = a + its 
norm is 

(1) N[x) = + 2}? + 4c^ - 6a6c = x ■ x' ■ x" 

with 

x' = a\ujh\pl\u?c\fl, 
X = a + a; 
and = a third root of one. 

Division in general leads to the corresponding field Q[v^]- Carrying 
out the rationalization of the denominator as in 

1 _ x' ■ x" _ - 26c + (2c^ -ah)^+ {b^ - ac)^ 

X ~ N{x) ~ N{x) 

gives the elements of Q[v^] in the form a + h\/2 + c\fi with a, 6, c 
rational fractions. 

3. The ambient vector space V 

As already mentioned, instead of (1, v^, v^) we use the algebraic basis 
(p, 1, 0") and elements of Z[v^] are expressed as 

w = x-p\y-Wz-(j 
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or changing the basis 

a + + c^fl = c ■ p + [a - 2c + h) ■ I + (h - c) ■ a. 

Now, let = M'^ be the 3-dimensional space endowed with the usual 
scalar product (a, b) and cross product a x b. Vectors can be written 
as triplets 

V = {v = (x, y, z);x,y,z e R} 
and define a linear mapping 

r] : Z[v^] ^ V 

by ri{x ■ p + y ■ 1 + z ■ a) = {x,y, z), the resulting image consisting of all 
vectors with integer entries, multiplication inherited from Z[\/2]. 

Taking into account 

x-p + y- l + z- a= {x + y — z) + {x + z)\/2 + x\/A 

we can define the norm function in the whole V: 

N{x, y, z) = {x + y — z)^ + 2{x + z)^ + 4x'^ — 6(x + y — z){x + z)x 

Multiplication with a will prove very important and we observe 

rj{a ■ w) = Sri{w) 

where 

1 
S= 10-3 
1 -3 

which, by the way, represents a hyperbolic toral authomorphism [5]. 
4. Interplay of different bases in V 

The term basis comes in several ways in mathematics. Our discussion 
needs it in two appearances: 

• as the basis of a number system [p in our case), 

• as the basis of a vector space (different bases of V here). 

The first usage figures in our paper [H], here we are concerned with 
the second one. 

We can think of the basis Bq = ((1, 0, 0), (0, 1, 0), (0, 0, 1)) as the canon- 
ical one. And we shall also write Bq = (s_i,So,Si) as we shall denote 
r]{a^) = Sj, and noting Ssj = Sj+i we have 

So = (0,1,0), 

si = (0,0,1), 

52 = (1,-3,-3), 

53 = (-3,10,6), 
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54 = (6,-21,-8), 

55 = (-8,30,3), 

56 = (3,-17,21), 

and we may also need the ones with negative indices 

s_i = (1,0,0), 

S-2 = (3,3,1), 

s_3 = (12,10,3), 

s_4 = (46,39,12), 

s_5 = (177, 150,46), 

s_6 = (681,577,177), 
as well as the inverse matrix 

^3 1 O" 



3 1 
1 



which we also meet in the Jacob i-Perron algorithm [T]. 



— (sj_i, Sj, Sj+i) for all integer 
j. These are good bases for our purposes as 



Further we define the series of bases Bi 

o 

Lemma 4.1. Elements ojl? have integer coefficients in each basis Bj. 



Proof. Let v be an element from 1? . If we multiply expansion 

V = asj_i + I3sj + 7Sj+i 
by integer element matrix 5*"-', we get vector with integer components 
S~^v = as_i + /3so + 7Si = (a, /3, 7). 

□ 

Besides the series of bases Bj we also define the conjugate series B* in 
the following manner. 

For start Sq = (1,0,0) and with the adjoint matrix 

'O 1 

^* = 1 
1 -3 -3 

we define vectors s*_|_]^ = S*s* for positive and negative indices. So we 
have 

s*_3 = (46,12,3), 

s*_2 = (12,3,1), 

s*-i = (3,1,0), 
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S5 = (1,0,0), 
St = (0,0,1), 

s* = (0,l,-3), 
sS = (1,-3,6). 

Remark 4.2. Comparing with the vectors Sj = {xj, yj, zj) we see that 

= (^j-i, = Zj-i, Zj). 

The bases series B* = (s*_i, s*, s*_^_j^) again providing integer coefficients 
for Zi[\/2]. The two bases series shall be useful in further computations. 

The action of linear transformation S is best understood in terms of 
its eigenvalues and eigenspaces. One eigenvalue is real, smaller then 1, 
and two are complex conjugate greater than 1 

Ai = a = -1 + ^, 

X^ = a' = -l + uj^=-l-f + 1^3^, 
A3 = a" = -1 + 00'^ = -1 - f - 
The eigenvectors being h and g ± ik where 

h= 1(^,2-2^+^,-^+^) = (0.209987,0.177926,0.05458), 
g = -^(_^,4+2^-^, ^-^) = (-0.104993,0.411037,-0.02729), 

k = ^{\/2, -2^-^, -^-^) = (0.181854, -0.592829, -0.410976). 
We can also compute the rotation angle 

9 = 71- arctan V- = 146.2° 
2 + ^ 

and after some computation we can express 

Sj = a^h + 2p2(gcos(j6') — ksin(j6')), 

and 

Sg = y/p{g cos — k sin 6*) , 
Sk = ^/p{g sin 6^ + k cos 9) . 

The norm ([T]) takes zero value on the union of the eigenplane P spanned 
by vectors g, k and the eigenline of h. Except for the origin, there is 
no rational point (x, y, z) of zero norm. 

The basic vectors Sj with increasing positive j are approaching the in- 
variant plane and for negative j being almost colinear to the eigenvector 
h. 

On the other hand we can construct the eigenbasis Be = (h, g, k) and 
the conjugate eigenbasis. 
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To make the conjugate eigenbasis, we compute the vector products 

h* = g X k, 
g* = h X k, 
k* = h X g 
and they constitute the conjugate eigenbasis 

h* = ^(1+^+^, 1, -1+^) = (-0.185104, -0.0481125, -0.0125055), 
g* = ^(_2+^+^, -2,2+^) = (-0.0407668,0.096225, -0.156843), 
k* = ^(i^ - 1^, 0, v^) = (-0.02729, 0, 0.104993). 
Later we shall need also the mixed product 

[h,g,k] = ^ = -M - -0.0481125. 
The plane P* of the vectors g*,k* is the invariant plane of S* . 

Lemma 4.3. The vectors h*,g* ± ik* are eigenvectors of the matrix 
S*. 

Proof. First we find the scalar products: 

(5*h*,g) = (g X k,^g) = (g X k, ^(gcos^ -ksin^)) = 0, 
(5*h*,k) = (g X k,5k) = (g X k, v^(gsin^ + kcos^)) = 0. 

So it is clear S*h* is orthogonal to both g and k and therefore colinear 
to h* itself, therefore an eigenvector with a real eigenvalue, the only 
one being cr, therefore S*h* = ah*. The other two cases demand a 
little more work to tell apart the two complex eigenvalues, of course 
unless we want to go into direct computation. □ 

We can also express the vectors s* in terms of the conjugate eigenbasis 

s* = -2V3^ [a^h* + pi (g* cos(ie - |) - MjO - |))) 

and infer a connection between the basis and conjugate basis via the 
matrix 

^3 1 O" 
1 
1 



namely 



and also connecting the eigenbasis with conjugate eigenbasis 
Th = -2VS\/2h*, 
rg=^(-g*-V3k*), 
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Tk= ^^{V3g* - k*). 
We shall also need some scalar and cross products of the basis vectors. 

Lemma 4.4. The scalar products of basis vectors are as follows: 
with 

{Sq, S_2) = 3, 
(So,S_i) = 1, 

{s*Q, So) = 0, 
(s*,si) = 0, 

(SS>S2) = 1, 

(s*,S3) = -3. 
Proof. Since 

we only need to read off the first component of Sj as Sq = (1, 0, 0). □ 

Lemma 4.5. For two consecutive basis vectors we have the cross prod- 
uct 

S-j X S_j+i = s* 
and if we jump by one index 

s_j_i X s_j+i = -s*_i + 3s*. 

Proof. Setting unknown coefficients a, (3, 7 

s-i X s_j+i = as*_i + /3s* + 7S*+i 
and taking scalar products in turn with s_j_i, s_j, s_j+i we get 

1 = a • 3 + /? • 1 + 7 • 

= a • 1 + ^ • + 7 • 

= a • + /3 • + 7 • 1 

and from here a = 7 = and (3 — 1. Likewise we prove the second 
formula. □ 
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J 


a 




8 


(20,-69,-33) 


5.34 


9 


(-9,27,20) 


1.68 


10 


(0,-7,-9) 


0.39 


11 


(-7,-9,0) 


0.28 


12 


(-30,-21,-7) 


0.65 


13 


(-111,-97,-30) 


1.89 



Table 1. Appropriate vector 



5. The shortest coefficient vector and convergents 

The vector (0,p, — g) = i]{p — qa) can as any vector be expanded in any 
basis Bj and with integer coefficients. 

In [11], using p basis, we expressed (0,p, —q) with a wider choice of vec- 
tors, but coefficients being hmited to 0,1,2 or 3: (0,p, —q) = YTj=k ^j^j- 
Whereas here, since each time only 3 vectors form the basis, we allow 
all integer coefficients. 

Should I be a convergent to a, we can control the size of the these 
coefficients, provided j has been chosen appropriately. 

Definition 5.1. Let | be a convergent to a and a = (01,02,03) the 
coordinates of ri{p — qa) = {0,p, —q) = aiSj_i -|- a2Sj + a^Sj+i in the 
basis Bj. This basis is called appropriate for the convergent, when the 

vector a' = p~4a is the shortest for some positive integer j. 

Example 5.2. From the table of convergents [TT] we take p = 1251, q = 
4813 that is just preceding the relatively big partial quotient bu = 14, 
so that |5| < ^ in the estimate p — qa = ^. 

Which j take, to make the vector a' shortest? Here are some results in 
the Table 1. 

And we see that the appropriate j and Bj to give the shortest a' is 
J = 11. 

Theorem 5.3. Let ^ be a convergent to a, Bj its appropriate basis. 
Then for its reduced coefficient vector we have 

|a'| < 2.01. 

Proof. We can write p = qa + ^ with |5| < 1 
(0,p, —q) = pso — qsi = p(h + 2g) — g(crh -|- 2y/p{gcos9 — ksin6')) 
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and when we rearrange the terms 

(0,p, -q) = -h + q^VsiVSg + k) + — g, 
q q 

we see what happens to either term under action of with grow- 
ing j. The first term grows exponentially with in the direction of 
the eigenvector h, the second decreases with p~i and rotates in the 
eigenplane, and the last term also decreases 



a^S-'{0,p,-q) = /7?-h + p-ig2^V3(gcos(j^ + ^) + ksin(j^+ ^)) 
(2) + (gcos(j^)+ksin(j^)). 

But if we put Kj — qp~^ we get 

ai = S-'{0,p,-q) = pi— h + pti^,-2^x/3(gcos(j^+ J)+ksin(j^+ J)) 

J\j D D 

2(5 

+ p"^ — (gcos(j^) + ksin(j^)). 

First we represent the vector a' = p"4a = a" + a'" with 

a" — + 2Kj v/2\/3(g cos aj + k sin aj), 

5 ■ 25 

a'" = — (gcosj6' + ksin j6'), 

where — j9 + ^. As for a" we write the square of its norm as a sum 
of three terms 



la'T = 



gcosccj -t-ksinccjl^ + g cos aj + k sin aj 



KJ 



= T1 + T2 + T3. 

The last term is estimated independently of Kj 

(3) iTal < 4 • 1 • ^V3max(h,gcosa + ksina) < 0.894896. 

Denote 

F{Kj, aj) = Ti + T2 = + 12\^K]\gcosaj + ksina^l^ 

j 

and to eliminate the dependence on a, and S we define another function 



J 

2 

2 



|h| 3 

G{K~) = —TT + 12v^i^f max |g cos a + ksin a 
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or inserting the numerical values 

G{Kj) < Q-Q'^^^^^^ + 12.95559953X| = H{Kj), 

j 

observing that F{Kj,aj) < H{Kj). The values of both functions de- 
pend only on the choice of j, the variable Kj assumes discrete values 
from a geometric series as Kj = qp~^. The function H{x) = ^ + hx^ 
as a function of continous variable a; > features just one minimum at 
Xo = ^fj^ with value H{xo) = but the discrete variable Kj shall 

almost certainly miss this minimum point. We shall further denote by 
x' the unique solution to the equation H[x') = H[x' p~* ). Indeed the 
equation H{x') = H{x'p^^) reads 

- + OX = ^ + ox p ^ 

x^p- 

and we can easily solve it 



with the value 

I 3 3 

E{x') = vab{p^ + 

However within the interval [x'p~^,x'] there is exactly one Kj and this 
defines also the choice of j. 

To determine j we have 

, _ 3 , 

X p 4 < Kj < X , 



, _3 _2i , 

X p t < qp t < X 



Taking logarithms 



3 3j 

\nx' — 7 In p < In (7 In p < In 

4 4 



and dividing by (— | Inp) we get 

, , , 4 In 4 In Q 4 In 

4 1-^^ >3-V^ — >-^^ ' 

6 Inp 6 Inp 6 Inp 

4(lng — Inx') 4(lng — Ina;') 

^ ^ ^1 > > ^1 ■ 

6 Inp 6 Inp 

4(lng — \nx') 



1 + 



3 Inp 



e N. 



Thus we have 

F{Kj,aj)<H{Kj)<H{x'). 
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Inserting numerical values gives Xq = 0.279205, x' = 0.462761, x'p 4 = 
0.168457, H{x') = 3.142064 and so F{K, a) < 3.142064, which together 
with the estimate ^ yields 

|a"| < V3. 142064 + 0.894896 < 2.009219. 
Let take e = 2.01 - 2.009219 = 0.000781. We approximate a'" 



I ///I -52 2 ■ 1 

I a \ < p 4 max 



5j 



; cos + k sin 9\ < p 4 



1.649395 

x'p~ 



< e 



from where we get j > 6 and the desired inequality follows for these j. 
If j < 5 we get from the inequality Q condition on q 



41ng 
31np 



> 



J 



4 In (7 4 In X 

> 

3 In p 3 In p 



and q has to be smaller than 73. There is only five convergents with 
such q and from the Table 2 we see, that computed |a'| satisfies our 
inequality. 



p/q 


|a'| 


1/3 


1.151 


1/4 


0.581 


6/23 


0.928 


7/27 


0.870 


13/50 


0.415 



Table 2. First five convergents 



Thus we can confirm 

|a'| < 2.01 

□ 

Remark 5.4. Over the first 10000 convergents we numerically find 
that |a'| < 1.753. The adjacent Figure [T] shows the statistics in dots, 
I a' I, for these convergents. 

Also the other way round, if coefficients are small enough, we are deal- 
ing with a convergent. The following theorem is however a rather coarse 
one. 

Theorem 5.5. Let a = (01,02,03) € be the coefficient vector for 
the basis Bj, j G N such that (a, s*) = and |a| < |p4. Then the 
resulting vector aiSj_i + 028^ + ass^+i = {0,p, —q) yields a continous 
fraction convergent -. 
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Figure 1. Length of vectors a' 

Proof. We multiply the equation ^ with the conjugate vectors h*, g*, k* 
to obtain 



(5|h*,a) 
(6|g*,a) 
(7r,a) 



p^-(h*,h) 



p-iq2^V3{g*, g) cosUe + ^) + p"^ -(g*, g) cosijO) 

o q 



26 



P 



'iq2\/2V3{k*,k) sm{j9 



■ p-5^(k*,k)sin(j^^) 



6 



Inserting the conditions of the theorem, we can estimate the quantitie 
— using (5): 



9 

(8) 



\5\ 3.|h*| 

q ^ 3M 



Using Q and Q we get 



q ■ 2^y3cos(je + -)M = pi(g*, a) - 2-Mcos{j9), 



q-2\/2V3sm{j9+-)M = -p2(k*,a) -2-Msin(j^). 



6 



We square and add up the last two equations to eliminate the sines 
and cosines 



q\2^V3Mf = p^((g*, a)2+(k*, a)2)+pi -M((k*, a) cos(j^^)-(g*, a) sm{je)) + ^M^ 

q 
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and estimate 



fi:i^uf < piMt^^ + i^(|g.| H. |k.|) H- p-^^i"-'^ 



2 



9 g MC, , , ,/ , g 

Using ([s]) and last inequality we can estimate 
|h 



5| < '" ' A/|g*p + |k*|2 + 4p~^|h*|(|g*| + |k*|) + 4p-3j|h*|2 



which is smaller than 0.48 for j > 2. For j = 1 condition |a| < |p4 < 
0.47 implies nonexistance of such integer vector. So we show \b\ < ^, 
that is enough for our conclusion. □ 

From the proof we can also infer the implications of an even smaller 
|a| on the \6\ and consecutively the next partial quotient B. 

Lemma 5.6. If in the above theorem |a| < ^p4,A < 1, we can esti- 
mate the next partial quotient B to the convergent | hy 



Proof. From the well known estimate [TO 

P 



1 

< 



we find < \5\ or 



< 



Bq^ 



But from the above proof, if instead of |, we put y, we have \5\ < 
■ 0.48 < ■ ^ and the estimate from lemma follows. □ 

Example 5.7. In our numerical experiment we found for j = 750, 
A = 0.03906, next B = 4941 = fegig (well known big partial quotient 
PTj ) with lemma suggesting B > 1308. 

5.1. The problem of the shortest lattice vector. Our case the 
lattice Aj = {a G Z'^,a ± s*}. Gauss reduction process mimics the 
euclidian algorithm. Let's have a basis zi,Z2, such that |zi| < |z2|. 
Choose k so that 

-^|zip < (Z2 - /i;Zi,Zi) < ^|ziP 

SO A; G Z is the nearest integer to -^r^^. 

o |Zl|^ 

Now, set the new Z2 := Z2 — kzi and compare: if |zi| < |z2| the process 
terminates, our shortest vector is zi, else we interchange zi -H- Z2 and 
start again. In some steps we get the shortest lattice vector [T4]. 
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Remark 5.8. In higher dimensions the so called LLL-algorithm |14j . 
[7] is similar to Gramm-Schmidt orthogonalization to generalize the 
Gauss process. 

Example 5.9. Here is how we carried out this process for j = 7, i.e. 

zi = s_6, Z2 = s_7 and we have the shortest vector a = (—7,1,0). 
Combining — 7s6 + IS7 = (0,59,-227) we read off the quotient ^ 
which does appear in the sequence of approximants (Table 3). 



n 


Zl 


Z2 


|zi|^ 


(Zl,Z2) 


k 


1 


(681,577,177) 


(2620,2220,681) 


828019 


3185697 


4 


2 


(-104,-88,-27) 


(681,577,177) 


19289 


-126379 


-7 


3 


(-47,-39,-12) 


(-104,-88,-27) 


3874 


8644 


2 


4 


(-10,-10,-3) 


(-47,-39,-12) 


209 


896 


4 


5 


(-7,1,0) 


(-10,-10,-3) 


50 


60 


1 


6 


(-3,-11,-3) 


(-7,1,0) 


139 


10 






Table 3. Shortest lattice vector 



We carried out the shortest vector algorithm for j = 2 until j = 1000. 
The resulting |a'| were, as shown in the dotted Figure |2| all below 
1. Only 21 of them did not result in continued fraction approximants 
(marked with squares). 



1.0 




Figure 2. Length of vectors a' resulting from the short- 
est vector algorithm 



6. X^-TEST OF DISTRIBUTION OF PARTIAL QUOTIENTS 

We applied x^-test to compare observed frequencies of partial quotients 
of with theoretical frequencies P{bn = k) = logj ^^r^- Using |T5] 
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we computed 75 000 partial quotients and divided them into R groups 
consisting of numbers 1, 2, 3, . . . , i? — 1 and of all numbers over R — \. 
Let Oi be the observed frequency of the i^^ group and Ei its expected 
frequency. The value of the test statistic is 



R 



{0^ - 

E, 



If the partial quotients the hypothesized distribution, X"^ has, approx- 
imately, a distribution with R—1 degrees of freedom. The resulting 
P-values for different degrees of freedom are shown in the Figure [3j 
Since all the P-values are above 0.05 we can not reject the hypothesis 
that the partial quotients of \/2 follow the distribution law of Kuzmin. 



1.0 r 



0.8 - 



0.6 - 



0.4 



0.2 




0.0 L 



df 



20 



40 



60 



100 



Figure 3. P- value for different degrees of freedom 
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